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Abstract 

The motion of a single charged particle in the space outside of a compact region of steady 
currents is investigated. The charged particle is assumed to produce negligible electromagnetic 
radiation, so that its energy is conserved. The source of the magnetic field is represented as a 
point multipole. After a general description, attention is focused on magnetic fields with axial 
symmetry. Lagrangian dynamical theory is utilized to identify constants of the motion as well as 
the equations of motion themselves. The qualitative method of Stormer is used to examine 
charged particle motion in axisymmetric multipole fields of all orders. Although the equations of 
motion generally have no analytical solutions and must be integrated numerically to produce a 
specific orbit, a topological examination of dynamics is possible, and can be used, a la Stormer, 
to completely describe the global aspects of the motion of a single charged particle in a space 
with an axisymmetric multipole magnetic field. 


1 . Introduction 


The problem of determining the motion of a single charged particle in a magnetic field is an old 
one but remains very challenging. The fundamental difficulty is that the equations describing the 
motion are generally nonlinear and no analytical solution is possible - computation is required to 
find specific trajectories. In the years 1907 to 1930, Carl Stormer and his small research group at 
the University of Oslo investigated the motion of charged particles in the geomagnetic dipole 
field, since this was important for understanding the behavior of solar particles and cosmic rays 
near the earth 1 . Although exact solutions to even this basic problem are not generally possible, 
Stormer found a method for describing global aspects of charged particle motion in a dipole field. 
His method was to look at qualitative solutions, whereby regions of space accessible to a charged 
particle could be mapped out in terms of certain constants of the motion. One constant is the 
energy (the charged particle is assumed to produce negligible electromagnetic radiation, so that its 
energy is conserved). The other invariant is related to rotational symmetry about the dipole axis 
and is called Stormer 's integral. Qualitative solutions which describe the partitioning of space 
into a collection of allowed and forbidden regions, with regard to particle motion, are themselves 
part of a larger, more modem discipline called topological dynamics 2 . (Well-known examples of 
spatial partitioning are the Van Allen radiation belts.) 

Here, the method of Stormer is extended and combined with other methods to investigate motion 
in axisymmetric magnetic fields other than a dipole field. In particular, point multipoles of all 
orders are considered. These sources are assumed to lie within a compact region of steady 
currents, so that any associated magnetic field tends to zero at large distances. At these distances, 
localized currents are well approximated by point multipoles, so that multipole expansions are 
the primary approximation far away from such sources. (However, if it is important to 
characterize charged particle motion near the actual currents that produce the axisymmetric field, 
then the magnetic fields - more specifically, the magnetic potentials - associated with current 
loops, rather than point multipoles, must be used; this topic will not be dealt with here.) 

Stormer's method for analyzing motion in axisymmetric magnetic fields essentially transforms 
the problem of finding specific three-dimensional trajectories into the problem of finding general 
two-dimensional topological solutions. 

An important aspect of the analysis presented here is the use of Lagrangian dynamical theory to 
identify constants of the motion as well as the equations of motion themselves. In particular, 
Stormer's integral for any axisymmetric multipole is easily found using Lagrangian methods. In 
order to do this, the relation between the magnetic scalar and vector potentials for point 
multipoles is developed. This result is necessary because Lagrangian dynamics explicitly uses 
the vector potential in its basic formulation, while it is the magnetic scalar potential that is 
usually expressed as a multipole expansion. Finding the “Stormer integral” for motion in the 
field of axisymmetric multipoles leads to the definition of an associated topological 
quasipotential. This, in turn, leads to an identification of the guiding center of the particle's 
motion as the absolute minimum of the quasipotential. 

The results developed are then applied to problems in astrophysics and in space propulsion 
system design. In astrophysics, qualitative solutions can be used to predict the shape of trapped 
particle regions in any multipole field (not just a dipole). In the design of space propulsion 
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systems that use strong, axisymmetric magnetic fields due to current loops, qualitative solutions 
can predict the ‘radiation belts’ that are bound to arise around spacecraft using these systems, 
while numerical solutions can produce the specific paths that charged particles take when leaving 
the spacecraft (for example, these particles can come from the edge of the exhaust plume of a 
plasma rocket engine or from interactions with the natural space environment). 


2. Lagrangian dynamics 

The basic equation of classical mechanics is Newton’s 2 nd Law: 



Here, p is the momentum of a particle of mass m. The most general form of p for a particle is the 
relativistic one: 



In the above, the velocity is v = dr/dt, where r is the position of the particle; v 2 = |v | 2 and c is the 
speed of light. 

A fundamental example of F appearing in (2-1) is the force due to an electric field E = -V<p 
derived from an electrostatic potential cp: 

(2-3) F = eE = -eV(p . 


The electric charge e can be positive or negative. The potential energy is V = ecp , so that 
equation (2-1) can be put into the form 


(2-4, a) ^ = 

dt dr 


b) f = V = e'^ie i -| T 

or dq i=i oq 


Here we introduced the Einstein summation convention: repeated indices imply summation. Note 
that, in the contravariant generalized components q' and basis vectors e', the superscript i is a 
coordinate index and not a power. There are also covariant quantities q, and e,. (The relation 
between the different bases is e' e y = 8) where 8) is the Kronecker delta. (For those unfamiliar 
with tensor analysis, please refer to the appropriate texts, for example, Sokolnikoff 3 or Lovelock 
and Rund 4 .) 

The momentum p can also be expressed as a gradient - the velocity- gradient of a quantity T: 
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(2-5) a)p = ^ b) T = -mc 2 7l-/v/ 2 /c 2 c)|- = e'^-. 

dv dv di? 

For v 2 /c 2 « 1, we have T + me 2 s {mv 2 , the non-relativistic kinetic energy. 


When the position vector is written as x, its components are specifically the Cartesian 
coordinates (x, y, z) = (x 1 , x 2 , x 3 ), with the associated constant orthonormal basis x,- , / = 1 , 2, 3, 

while the components of the vector q are the generalized coordinates, q‘, i = 1, 2, 3, and the 
associated independent basis consists of the vectors e„ / = 1, 2, 3. (The set of available q k 
includes the possibility that q k = x k .) The position vector r resides in position space, while the 
velocity vector v = dr/dt exists in velocity space. The components of q are q l = dq' j 'dt (in 
general, q ^ v ). The space consisting of vectors q is the configuration space for a single 
particle; the six-dimensional space (q, q ) is often called phase space. (In the case of n distinct 
particles, the full configuration space has 3" dimensions and the full phase space has 6" 
dimensions.) 

In phase space, derivatives can be taken with respect to both position and velocity. Consider the 
following derivatives: 


( 2 - 6 ) 


a) 


d£ 
dc . 'h 


d ( dx i .j 


\ 


/ 


dx‘ 

dq k 


d dx l _ av .j _ dx‘ 
dt dq k ~ dq j dq k 9 ~dq k ' 


These quantities are used in the coordinate transformations that follow. 

We take Newton’s 2 nd Law (2-4a), write it in Cartesian coordinates and transform to generalized 
coordinates, using the Jacobian dyadic (i.e., 2 nd rank tensor) J : 


(2-7) a) J — = -J — 
dt dx 


, s „ dx dx' k „ 

b) J = aTv e X ' 


Using (2-4b) and (2- 7b), the left side of equation (2-7a) becomes 
( 2 - 8 ) 


d 

r dT ' 

a a.v' d 

' ar ' 

~e k 

d 

( dx‘ dT " 

dT d dx‘ 

dt 


_e dq k dt 



dt 

dq k dx J 

di‘ dt dq k 


Then, using (2-6b) on the right side of (2-8) gives 
(2-9) J-- 


d 

f dT ' 

k 

d 

' dT ' 

dT 

dt 


= e 

dt 

kJ 

dq k _ 


The right side of (2- 7a) is simply 
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( 2 - 10 ) 


dy__dx k dV 

dx dq dx dq k 

Since the e k are independent, equation (2-4) takes the form 


( 2 - 11 ) 


d(dT' 

~di\W y 


dT 

V 


dV 

w 


Finally, in the case where the potential is V = e(p, where V is only a function of physical space r 

(i.e., only of the q' ), we can define a quantity L = T- V, called the “Lagrangian function” or just 
th e Lagrangian. In generalized coordinates and velocities, equation (2-1), then becomes 


( 2 - 12 ) 


d ( dL > 
*{dq k , 


dL 

dq k 


= 0 . 


This is the Euler -Lagrange equation*'*' 5 . In a Cartesian system q‘ e {x,y,zj and q‘ € {x,y,zj, 

while in a spherical polar coordinate system q‘ e {r, 0,<J)} and q 1 e {r,9,(j)} . In fact, the q‘ need 

not be orthogonal coordinates and the e, need not be orthonormal, just independent. Furthermore, 
equation (2-12) remains valid even if force is dependent on velocity, as well as position. 


The principle of least action 

Equation (2-12) can be shown to result from requiring that a certain integral, called the action, be 
a minimum with respect to variation of the generalized coordinates q‘ and velocities q' about the 
trajectory of a particle (allowing no variation at the beginning and ending times t\ and ti) The 
action S is (where q and q are shorthand for q‘ and q‘ , z = 1,2,3, respectively): 


(2-13) S = ‘]L{q,q)dt. 
1 1 


The development of the Euler- Lagrange equation from the principle of least action proceeds as 
follows 3 ' 4 ’ 5 . 


The principle of least action requires that the variation of the action S with respect to any 
independent variables be zero to first order. Let us vary the generalized coordinates and velocities 
within the action as given in (2-13), but not at the end points t\ and ty. 


(2-14) 8S = ‘}hL(q,q)dt = ){ |^-8 q‘ +-^-8 q‘ 


\ 


dq‘ 


dq 1 


dt = 0 . 
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Variation of the generalized velocity satisfies 


(2-15) dq‘ = d = 

dt dt 


Upon integrating (2-14) by parts and using (2-15), we have 


(2-16) 85 = ) 


^dL • dL ddq' 
-dq + — - 

Z n • l 


dq‘ 


dq dt 


1 , 

' 2 '}( 

dt = — -dq 

-) 

J dq 1 



d dL dL 


dtdq dq' 


dq' dt . 


Since there is no variation at the endpoints, we have dq l {t\) = 8 q‘{ti) = 0> so that the first term on the 
right side of equation (2-16) is zero. Finally, since dq 1 is small but otherwise arbitrary for t\ < t <tj, 
we see that (2-12) is required for dS to be zero, to first order. 

Lagrangian for a charged particle 

In electrodynamics 6,7 , we proceed as follows. First, in the case when F is due to an electrostatic 
potential (p, as in (2-3), the Lagrangian is 


(2-17) L = -mc 2 J\- / v / 2 / c 2 - 


e(p. 


The potential cp does not include a contribution from the particle itself, but is due to other 
charges; if there are no others, then the potential is set to zero (or some other constant). 

In general, a magnetic field B, as well as an electric field E, will exist. Then, in addition to the 
scalar potential cp, a vector potential A will be needed: 

1 d\ 

(2-18) a) B = V x A b) E = V(p. 

c dt 

Note that the divergence of the first equation and the curl of the second equation above produce 
two of Maxwell’s equations: 

(2-19) a) V • B = 0 b) VxE . 

c dt 

% 

For completeness, and for later use, the other two Maxwell’s equations are 


1 A — 

(2-20) a) V • E = 47tp b)VxB=-— + — j. 

c dt c 
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Equations (2-19) are known as the homogeneous Maxwell equations, while those in (2-20) are 
called the inhomogeneous Maxwell equations because they have as sources the charge density p 
and the current density j. For a single particle, p = ^5(x-x (1) ) and j = gVfi^x-xp)), where S(x) is 
the 3-dimensional Dirac delta function, x^, is the particle’s actual position and v (1) its velocity at 
time t. (In the present work, we will be interested, eventually, only in one restricted case: the 
magnetic field in a source- free region, where, except for the single charge itself, the current 
density j and charge density p are both zero.) 

In the many-particle case, the Lagrangian for a collection of« charged particles takes the form 
(the summation here is over particles rather than coordinate indices): 

( 2 - 21 ) 


L = ~ I 

i=i 


'{ m (i)C 2 ijl-vji/c 2 + e {0 (p(x (0 ) - — A(X (0 ) V (0 1+ — J \J E(x') | 2 - 1 B(x')/ 2 ] d V . 

1 ^ C j 871 


The potentials <p(X(,-)) and A(X(,)) do not include contributions from the z-th particle. The last term 
is the Lagrangian for the electromagnetic fields E and B; it is a constant under variation with 
respect to the spatial variable x,- because the integration is over all space. However, variation of S 
in (2-13), with L given by (2-21), with respect to 9 , the components of A, and all their 
derivatives (which serve the same roles for a continuous field as the q' and q do for a point 
particle), yields the electromagnetic field equations (2-19) and (2-20). This is a result of classical 
field theory 6 , which, although highly interesting and important, is outside the scope of this work. 


3. Lagrangian of a particle in a steady magnetic field with axial symmetry 

The Lagrangian for a single particle in a steady magnetic field 

In the present work, the primary concern is with only one particle in an externally imposed, time- 
independent magnetic field (i.e., no electric field E is acting on the particle). The Lagrangian is 
given by (2-21) for one particle (n = 1). We also set cp = 0 and omit the last term (which is 
merely a constant with respect to the particle’s generalized coordinates and velocities, q, and q t ). 

In this case, (2-21) becomes the effective Lagrangian for a single charged particle in a steady 
external magnetic field (with x = x (1) and v = v ([) ): 


(3-1) Z, = -mc 2 i/l-v 2 /c 2 +eX (\)\/c. 

The constant magnetic field B = Vx A may be found by using (2-18) along with 9 = 0 and the 
fact that A is not a function of the time t. 

We also assume that any radiation fields that the particle may produce are negligible. Otherwise 
we would need to include the last term in (2-21), as A would no longer represent a steady (i.e., 
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time- independent) magnetic field. A single particle in a steady magnetic field has an important 
constant of the motion: its energy. This is shown by developing a general result. 


The first constant of the motion 

To begin, we sum the product of q t with the expression found in equation (2-12): 


(3-2) 


.id 
q — 
dt 


(ail 

, i dL 

d 

f.tdL) 


f..,- dL dL \ 

dq' 

’ J 

q dq 1 

dt 

q dq 1 
\ ^ ) 


q +q ^ < 
O? 3? J 


= 0. 


Above, we have differentiated by parts (and again used the summation convention). The quantity 
in the last parenthesis on the right of (3-2) is the total time derivative of L : 


„ „ dL dL . f dL 
(3-3) — = q ‘— + q' — 

dt dq dq 


If L were an explicit function of the time, then it would be necessary to add the term dL /dt to 
the right- hand- side of (3-3). 

Next, using (3-3), equation (3-2) becomes 


(3 - 4) n 


• i dL 
q — : — L 

dq' 


\ 


= 0 . 


The conserved quantity in parenthesis is the energy H : 

(3-5) H = q'^--L. 

dq 


Thus, once the form of the Lagrangian (3-1) is specified, the energy H immediately follows. 

Equation (3-5) is a general result that depends only on the fact that L is not an explicit function 
of time. [Although we will only consider cases with only a magnetic field present, an electric 
field could also be included. As long as both electric and magnetic field are steady, then the 
energy H, as determined from (3-5), still remains constant during a particle’s motion.] 

Axially symmetric magnetic fields 

Now, consider the case where the magnetic field is axially symmetric or axisymmetric, that is, 
invariant upon rotation about some axis. In this case, two appropriate coordinate systems are the 
spherical polar coordinate system ( q t e {r, 0, <J>} ) and the circular cylindrical coordinate system 
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(<7, e {p,<|), z}). Specifically, it will be assumed that the magnetic field does not change when ([) 
varies, i.e., B is not a function of <j> (or t). 

In what follows, we will see the usefulness of the Lagrangian method, as embodied in (2- 1 4), 
when L is not a function of the coordinate 0. (Here, the absent coordinate is <|>, but in many other 
physical situations, it may be some other generalized coordinate <?,.) The importance of a missing 
coordinate in the Lagrangian is that this leads to a second constant of the motion (as will be seen 
in the next section). 


4. Lagrangian with axial symmetry in spherical polar coordinates 

In a spherical polar coordinate system, the coordinates are <7, e {r,0,(j)}and an axially symmetric 
magnetic field has the form: 

(4-1) B = B r (r,Q)r + B e (r, 0)0 <=> A = J(r,G)+. 

In the above equation, r, 0, ^ are the unit vectors in the r, 0, (f) directions, respectively. 


Using A as defined in (4-1), the Lagrangian (3-1) becomes 
(4-2) L = — me 2 -^1 - v 2 /c 2 + eA(r,G)v (Sf jc . 

Here, v 2 = v r 2 + ve 2 4-v,), 2 , with v r = r , v e = rQ and s r sin 0<j) . Since L does not contain <|> (it is 

said to be “cyclic in ([)”), then the <|) component of (2-14) can be immediately integrated to 
produce a constant quantity: 


(4-3) 


d(d£ 
dt [ 5<i> j 


= 0 


— = A = constant . 

a<j) ♦ 


We can now use (3-5), (4-2) and (4-3) to produce the explicit form of the constants of the motion 
for the case of an axially symmetric magnetic field. First, using equations (3-5) and (4-2), as well 

as = rsin0<j) , gives an expression for the constant energy H: 


(4-4) 


H = 


me 


-i/l - v 2 / 1 


Second, placing (4-2) into (4-3) yields the constant conjugate momentum P$. 
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(4-5) = r sinG 

9<t> 


mrsim 




+-A{r,d) 


= P, 


Equation (4-4) can then be used with (4-5) to write in the form 
rsin 0 


(4-6) 


n = ' 


H 

— rsin0<}> + e/l(r,0) 

c 


Thus, in the case of an axisymmetric, steady magnetic field, there are two constants of the 
motion: the energy H, as defined by (3-5), and P§, known as Stormer’s integral. 


5. Equations of motion in spherical polar coordinates 

In this section, exact and qualitative equations of motion will be discussed. 

Equations for a single trajectory 

The equations of motion of a charged particle in an axially symmetric magnetic field, specified 
by A(r,Q), are found by putting the Lagrangian (4-2) into the Euler- Lagrange equations (2-12), 
with q t € {r,0, <)>} ; the equation for (j) comes from (4-6). The results (remember that H and P§ are 
constants) are: 


(5-1) 


. dr ecdirA ) . v 2 -r 2 dr 

a) ~T = ~77~\ — -sm0<|>+ b ) ~ = r 

dt H dr r dt 


(5-2) 


a) 

dt 


ec d (sin 6 A) 

77 ao 


+ /-sin0cos0<j) 


b) — = 9 
dt 


(5-3) i = 

dt r sin0 


c 

L 

7/v 

r sin0 


-eA(r,6) 


The velocity magnitude v is introduced into (5-3) for convenience later. Some mathematical 
manipulation 6 can show that (5-1), (5-2), and (5-3) are equivalent to p = evxB/c. 


Equations (5-1) - (5-3) comprise five highly nonlinear equations that serve to determine the 
exact trajectory of a charged particle in an arbitrary axisymmetric magnetic field B(r,0) with a 

vector potential A = ^(r,0)4, once a set of initial conditions are given. One constant of the 
motion, Stormer’s integral, is used to produce equation (5-3), while the other constant of the 
motion, the energy H, can be used to reduce the number of independent equations from five to 
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four. It is well known that such a set of equations describes chaotic motion 1 1,12 , i.e., if the initial 
conditions are changed only slightly, then the resulting trajectory diverges quickly from the one 
flowing out of the original initial conditions. In chaotic motion, no exact analytical solution is 
generally possible if A is nontrivial - a numerical solution is required. However, a qualitative 
solution is possible, i.e., one that allows a mapping of the regions accessible and inaccessible to 
particle motion. We now turn to this topic. 

Qualitative solutions 

In the following, we use well-established (though not necessarily widely-known) results 1,11 that 
were developed by Stormer in his studies of the polar aurorae. These results allow a qualitative 
solution of the equations of the highly non-linear motion of a charged particle. 

Since H= mc 2 (l-v 2 /c 2 )' l/2 is constant, so is v 2 . In a spherical polar coordinate system, v 2 has the 
form 


(5-4) v 2 = v 2 + Vq + = r 2 + r 2 0 2 +r 2 sin 2 0(j) 2 . 

Now, if we divide the expression in (5-4) by v 2 , and use ds = vdt ( ds is an element of the path 
length along the particle’s trajectory), along with (5-3), we get 


(dr) 

2 

( 

— 

+ 

y 

ds 


ds 

V ) 


\ J 


This can be thought of as the total energy equation for a surrogate particle 9 ' 12 - one of mass 2, 
moving in a two-dimensional potential energy field h 2 . 

Note that since v is a constant, so is p = [p|, where p is given by equation (2-2). Using pc 2 = vH 
and (5- 3b), we see that h can be expressed as 

(5 . 6) 

prsin 0 pc 


We will call the quantity h 2 the pseudo-potential. We must have 0 < h 2 <1 during a particle’s 
motion, as (5-5) shows, so the motion of a charged particle in an axisymmetric magnetic field is 

restricted to those areas of the r-0 plane for which — 1 < h < 1 ; these are allowed or accessible 
regions for particle motion. Regions for which h 2 > 1 are forbidden or inaccessible regions This 
fact will allow a qualitative global analysis of charged particle motion once a functional form for 
A(r,Q) is chosen, a topic which will be addressed next. 
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6. Multipole Fields 


In the case of a steady magnetic field in a space r > r 0 where there are no sources and no electric 
field, equations (2- 19a) and (2-20b) yield the following field equations for B: 

(6-1) a) V B = 0 b) V x B = 0 . 

Thus, B can be expressed as either 

(6-2) a) B = V x A or b) B = -V<t> . 

Equation (6-2a) was already seen in (2- 1 8a), where A was identified as the magnetic vector 
potential. In (6-2b), O is termed the magnetic scalar potential. If we place B = -VO into (2- 
19a), i.e., V • B = 0, we obtain Laplace’s equation: 

(6-3) V 2 O=0. 

Requiring that O — » 0 as r — > in a spherical polar coordinate system leads to the following 

well-known general solution 7,8 : 

(6-4) <t> = £ £ -^-iUe.4>) . 

n=\ m=-n V 

Here, the c nm are complex constants satisfying c* nm = (- (the * indicates complex 

conjugation), and the Y„ m are the well-known spherical harmonics : 


(6-5) 7„ m (0, <(>) = . 


| 2 « + 1 («-/«)! 

4n (n + m)! 


p,:( 


COS' 




The associated Legendre functions P"‘ satisfy 


m jm 


(6 - 6) k w = t 2 - 1 )" = (- 1) m (i - x 2 ) ^ — p n w . 


The P n (x) = P„(x) are the Legendre polynomials. Furthermore, we have 

(6-7) a) P„-"(xH-irj^^P„ w (x) b) 7„,_ m (0,())) = (-l) m 7; m (0,(()). 

\n + m)\ 


In Pf , n is the degree and m is the order the of the associated Legendre function. 
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The vector potential equivalent to a scalar potential 


In (6-2), two possibilities are given: B = Vx A or B = -VO . Here, we will find the vector 
potential A,,,,, equivalent to the scalar potential 0,„„, where 0,„„ is one of the terms in the 
spherical harmonic expansion (6-4) of O : 


( 6 - 8 ) 


a) O = XO, 


b) <*>„„, = 


f 

J tli 


cose)e 


imt j) 


First, we will assume that A has the fonn 

(6-9) a)A=ZA„ b) A„=%[F-(e)e+o;(e)tle'”». 


Using (6-8) and (6-9), along with the requirement VxA nm = -VO /im , yields 


(6-10) ^-(smQG:)-imF n m = (« + 1) sin BP n m 

u u 


(6-11) a) „GZ=-1S- b) nF,:=~P: 

dB srnG 


Here, Z 5 ”' = Z 3 "’ (cosG) . Placing F™ and as given in (6-11) into equation (6-10) yields 

( 


( 6 - 12 ) 


1 d 


sinG dB 


sin 6 


dPr ^ 


d G 


+ 


n(n + 1) - — 


m 


sin 2 G 


o. 


V J L -» 

This is the differential equation satisfied by the associated Legendre function P,"'(cosG) . Thus, 
the three equations in (6-10) and (6-11) are consistent. 

Therefore, if the magnetic scalar potential O contains a term O n „, as given in (6-8), the magnetic 
vector A contains a term A,,,,, given by 


(6-13) 


\ f nm 

A nm 


nr 


n + 1 


imP„ m { cosG) a dP n m { cosG), 

® — f 


sin G 


dB 


Jmty 


Another way of writing this is 


(6-14) 


A = 

nm 

n 


e d 

sin 9 3(j) 


♦- 


= --rxVO 
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The vector potential, in terms of a given scalar potential, can also be expressed as 
(6-15) a) A„ m =Vx0„ ffl b) 0 nm =irO„ m . 

The vector potential A nm clearly satisfies V • A nm = 0 . 

Axisymmetric multipoles 

Now that we have the 3-dimensional form of A, we will reduce this to the form of A in (4-1), 
i.e., Aq = 0 ,Aq = A(r,Q), which is appropriate for a magnetic field with axial symmetry. In this 
case, we use (6-13) with m = 0, A„ = A, l0 , and a n =f n o- 


(6-16) a) A„ =A(r,e)+ b) A(r,6) = - % dP »( C ° sd [ 


nr 


The associated Legendre functions satisfy (among others) the following recursion relation: 


(6- 1 7) — P’" 1 (cos0) = P"‘ (cosG) + (m- l)cot 0F"' 1 (cos0) . 

dQ 

(6- 1 8) P™ +2 (cos0) + 2{m + l)cot0 P™ +1 (cos0) + (n - m){n + m - 1 )P” (cos0) = 0 . 

If we set m = 1, we can use (6-17) to express A (r,B) in (6- 16b) as 

(6-19) A(r,0) = — -~n+r^n (cos0). 

nr 

This is the n th order multipole component of an axisymmetric magnetic vector potential: n = 1 is 
a dipole, n = 2 is a quadrupole, and n = 3 is an octupole, etc. 


7. Lagrangian with axial symmetry in circular cylindrical coordinates 

In a circular cylindrical coordinate system the coordinates are q' e {p,0, z} and an axially 
symmetric magnetic field has the form: 

(7-1) B = (p, z)n + 5 Z (p, z)z <=> A = A(p,z)+. 

In the above equation, n, t z are the unit vectors in the p, p, z directions, respectively. 

Using (7-1), the Lagrangian (3-1) becomes (with v 2 = v, 2 + + v\ ; v p = p , = p<(), v z = z) 
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(7-2) 


L = -mc 2 ^j 1 - v 2 /c 2 +e^(p,z)v ()) /c. 


Again, since L does not contain 0, then the <J) component of the Euler- Lagrange equations (2-12) 
leads to 


(M > n 






dL 

= 0 => — = r fh = constant . 

* 


Putting (7-2) into (7-3) again yields Stormer’s integral, this time in cylindrical coordinates: 
(7-4) 


A =- 


H 

— pty + eA(p, z ) 
c 


The energy H of the charged particle is defined in terms of L by (3-5). After, putting (7-2) into 
the general expresson (3-5), we see that it takes the same form as (4-5): 


(7-5) 


H = 


me 


Vl — v 2 /c J ' 


Next, we will discuss the equations of motion in circular cylindrical coordinates, as it often 
useful to phrase the solution in terms of a cylindrical rather than a spherical polar coordinate 
system. 


8. Equations of motion in circular cylindrical coordinates 

Equations for a single trajectory 

The equations of motion of a charged particle in an axially symmetric magnetic field, specified 
by 4(p, z), are found by putting the Lagrangian (7-2) into the Euler- Lagrange equations (2-14), 
with q i € {p,(f), z} ; the equation for <j> is determined from (7-4). The results (again, remember 
that H and are constants) are: 


(8-1) 

. do ecdipA )■ 12 

dt H 3p 

dt 

(8-2) 

. dz ec d A 

a) — = — P— P 

dt H dz 

u\ dz 
b) — = z 

dt 
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(8-3) 


d§ _ ^ _ v h 
dt p 


b) h= 

vH 


cP* 

eA(p,z) 

P 


Again, some mathematical manipulation 6 will show that these are equivalent to p = evxB/c. 


Equations (8-1) - (8-3) are equivalent to (5-1) - (5-3) and also comprise five highly nonlinear 
equations, that serve to determine the exact trajectory of a charged particle, in an arbitrary 

axisymmetric magnetic field B(p, z) with a vector potential A = A(p,z )^, once a set of initial 
conditions are given. One constant of the motion, Stormer’s integral (7-4), is used to produce 
equation (8-3), while the remaining constant of the motion, the energy H as given in (7-5), can be 
again used to reduce the number of independent equations from five to four. These equations 
again define chaotic motion, in which no analytical solution is generally possible if A is non- 
trivial, leading to the necessity of numerical solution. 

Qualitative solutions 

Once again, a qualitative solution is possible, i.e., one that allows a mapping of the regions 
accessible and inaccessible to particle motion. At large distances from the assembly of current 
loops, the magnetic field and vector potertial can be represented in terms of multipoles and the 
discussion presented in Section 5 is again pertinent. At closer distances, the vector potential 
appropriate to a set of current loops is required (and will be provided in the next section). 

As before, v 2 is constant. In circular cylindrical coordinates, v 2 has the form 
(8-4) v 2 = v 2 + v 2 + v 2 = p 2 +p 2 (j> 2 + i 2 . 

Dividing the expression in (8-4) by v 2 and using (8- 3b), along with ds = vdt ( ds is again an 
element of the path length along the particle’s trajectory), gives 


( 


(8-5) 


dp ) 2 fdz^ 2 
+ 


V 


ds 


ds 


+ /i =1 


/ 


As with (5-2), this can be thought of as the total energy equation for a surrogate particle of mass 
2. Here, the two-dimensional pseudo-potential h 2 is a function of p and z. 

We have already seen that if v 2 is constant, so is p = |p|. Using pc 2 = vH, along with equations (8- 
3b), we see that the function h can be expressed in cylindrical coordinates as 


( 8 - 6 ) 


, d ( )> 

h =p — 

ds 


A, eAjflz) 
P P PC 


As before, 0 < h 2 <1 during a particle’s motion, so that the motion of a charged particle in an 
axisymmetric magnetic field is restricted to those areas of the p-z plane for which — 1 < h < 1 
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(these are regions where particle motion is allowed). A particle is forbidden to travel in regions 
in which h 2 > 1, while h 2 = 1 defines the boundary between the two regions. As in the case of 
spherical polar coordinates, this will enable us to perform a qualitative global analysis of charged 
particle motion once a functional form for A(p, z) is chosen. 


9. Global structure of the pseudo-potential 


The pseudo-potential h 2 will be expressed either in spherical polar or cylindrical coordinates, 
where h is given by (5-6) for spherical polar and by (8.6) for cylindrical coordinates. Since the 
pseudo-potential is the square of h, some general comments can be made on its global structure. 
Let (q\ = r, <72 = 0 or q\ = p, <72 = z); then the extrema of h 2 are found by requiring 

(9-1) — = 2/i— =0, i=l,2. 

a<7, dg, 

Obviously, h = 0 is an absolute minimum of h 2 . There is also the case where dh/dq; = 0. 

The nature of the extremum is found by considering the second derivatives of h 2 : (again, we use 
<71 =r, <72 = 6 or 91 =p, q 2 =z) 


(9-2) 


d 2 h 2 J. d 2 h aA dh 

= 2 h + - 


dgfaj 


V 


dqfiq - a<?, dq , 


i = 1,2; j = 1,2 . 


If h 2 is expanded in a Taylor series about the point q t = q ip + dq t and (9-1) holds, we have, to 
second order, 


(9-3) 


h 2 -h 2 = 


d 2 h 2 


dq 2 


2,2 


rh 

“<7i +2 — 

oq { oq 2 


d 2 h 2 


dq \dq 2 + _ 2 

a^ 2 


dq l 


If h 0 = 0, then (9-3) becomes 
(9-4) 


dh 

J dh 

dq 1 + _ 

) 

dq 2 

dq { 

O 

CL 

N) 

0 > 


h l = 2 

Using (9-4), we see that the equation of the line along which h = 0 is 
dq { _ dh 


(9-5) 


dq 2 dq 2 


( dh 


o[ a <?i 

O 
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A low energy particle will be confined close to this curvilinear minimum of the pseudo-potential. 
This is the so-called guiding center of the particle’s motion; more detail will be provided 
presently for the point multipole and current loop cases. 

Using (9-1), it follows that in the case when dh/dqi = 0, the discriminant D of the quadratic form 
on the right side of (9-3) is, 


(9-6) 


D = \6hl 






d 2 h 

d<lid<h\ 0 




J 


d 2 h 





\ 


J 


If D > 0, the point at which dh/dqi = 0 is a saddle-point, since the quadratic fonn on the right side 
of (9-3) will be divided into regions of positive and negative values. On the other hand D < 0 
signifies a local maximum or minimum, and D = 0 denotes a level line, similar to (9-4). 

Global structure - point multipoles 

Combining equation (8-6) with (6-19) yields h for the point multipole of degree n: 


(9-7) 


h =• 


-+- 


/?rsin9 pc nr 


Al+1 


K (cose) . 


We now define the Stormer radius rs and constant y, and the dimensionless variable R\ 


(9-8) a )r s n+ '=^L b) Y = —^ c)R = — . 

npc (n + l)pr s r s 

Using these (9-7) can be written in a clearly non-dimensional form: 

(9 -9) h = l£±])l + ^(£2!9). 

R sinG R nM 

Thus, the topology of the pseudo-potential h 2 is controlled by the dimensionless parameter y. 

The absolute minimum of h 2 occurs when h = 0. Although h 0 as R , using (9 -9) gives the 

relation between finite R and 6 for an absolute minimum: 

(9-10) a) (n + l)y/?p = -sinG (cosG) b) y/^cosG) < 0. 


Let us define L 0 "= (n + l)|y| and assume that (9- 10b) holds; then (9- 10a) takes the form 
(9-11) R” = To sinG |p,!(cose)| 
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Using the table in Appendix A, the case n = 1 becomes Ro = Z, o sin 2 0. As is well known, this is the 
equation for a field line in a dipole magnetic field - the guiding center of a particle’s motion in a 
dipole field. Only in the case n = 1 is it possible that (9- 10b) holds (with y > 0) over the entire 
range of angles: 0 < 0 < n. For those cases with n > 1, there are intervals between 0 and 7t in 
which 0 takes values such that the right side of (9- 10b) does not hold for either y > 0 or y < 0; in 
these intervals, no absolute minimum exists for finite values of R. 


Using (9-9) and (6-17), the derivatives of h are readily found to be 


(9-12) 


dh _ n + 1 _y_ + P^i cos0) 
dR R 2 sin0 R " 


(9- 13) ~ = - ( ”+ 1 >r^ te + _l_[^ (cos9) + cote Pl (cose)]. 


7?sin0 


R 


At extrema, the conditions dh/dR = dh/dB = 0 hold; these conditions define the coordinates R p 
and B p of the extremal points. Using equations (9-12), (9-13), and (6-18), we have 


(9-14) 



sinQ^ Pj (cosQp) 
Y 


(9-15) cos 0, P x n (cos B p )-{n + l)sin B p P n (cos 0„) = P„' +] (cos 0,) = 0 . 

[In (9-15), #8.733-3 of Ref. 12 has been used.] Solving these equations yields R p and B p ; since 
(9-15) contains a factor that is an « lh order polynomial in cosB p , there will be, in general, n 
different sets of real values for R p and 0 P with 0 < 0 P < 7 t . If a? is odd, B p = 7t/2 is a solution; the 
other n - 1 solutions come in ( n - l)/2 pairs Q p = a p and Q p = n - a p , with 0 < a p < n/2. If n is 
even, then there are nil pairs a p and 7t - a p , again with 0 < a p < 7t/2. For n > 1 , different values of 
B p will give the right side of (9-14) a positive value or negative value, so that for a fixed sign of 
y, not all solutions of equation (9- 1 5) lead to the location of an extremum, but only those for 
which the right side of (9-14) is positive. 

Placing (9-14) into (9-9) yields 


(9-16) 


n y 

R p sxrB p 


At this point, we can take the second derivatives of h with respect to R and 0, and evaluate them 
at R p and B p . Differentiating (9-3) and (9-13) results in 
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(9-17) 


d 2 h 

dR 2 


n{n + l)y 
RlsmQ p 


(9-18) 


d 2 h 


dQ 1 


= —3 [l + (/t + l)sin 2 0 ] 

P R P Sm e p 


(9-19) 


d 2 h 


dRdQ 


n(n + l)ycot0 / 
R 2 p sin0 p 


It is evident that the 2 nd derivatives in (9- 1 7) and (9- 1 8) are negative definite and positive 
definite, respectively, if y > 0, and positive definite and negative definite, respectively, if y < 0, 
and both zero if y = 0. The discriminant D, as given by (9-6), is thus positive if |y| > 0, and it has 
value D = 0 only if y = 0. Using (9- 1 6), (9- 1 7), (9- 1 8), and (9- 1 9), along with (9-6), gives 


(9-20) 


D = 16 


2 2 

n y 


( n + 1)(« + 2) . 


Since D > 0 for |y| > 0, R p and 0 P are the locations of saddle points, if the right side of (9-14) is 
positive. 


Setting n = 1, ..., 5 in (9-15) and using the (cos0 p ) from Appendix A, leads to the results in 
Table I (for -1 < cos0 p < 1): 

Table I. Values of cos0 p of possible saddle points for n = 1 to 5. 


n - 1 

COS0/, = 0 

n = 2 

cos0 p = ± a/5 

n = 3 

cos Q p = 0, ±a/3/7 

n = 4 

cos0 ^ = ± ^ 


n = 5 

COS0 P = 0, ± y 

's± 2 jTn 

11 


In those cases that have n > 5, equations higher than quadratic must be solved. For n = 6 and 7, 
the solution of a cubic equation is needed, and for n = 8 and 9, a quartic equation must be solved 
(although these will lead to closed form solutions). For n > 9, a numerical solution is generally 
required. Here, we stop at n = 5, as this is sufficient for our purpose. 
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In order that a particle be able to actually pass thru a saddle point, we must have the pseudo- 
potential at the saddle point satisfy \h p \ < 1. Setting \h p \ = 1 in (9-16), using (9-14) and the values 
of 0p for n = 1 to 5 in Table I, allows the critical value y c to be determined: 


(9-21) a) |y c | = sin 0 /^‘(cosG^) 


rt +1 pi rt + 1 


b) Y c P„ l (cos0 p )<O. 


Once y c is found, equation (9- 14) can be used to find R p . Using this procedure for n = 1 to 5, 
gives the results in Table II. The meaning of y c is the following: If |y| | Yc|> the particle can no 
longer pass through the associated saddle point. 


Table II. Saddle point positions and critical values y c . 


n 

R P 

Q p (degrees) 

Yc 

i 

1.0000 

90.0000 

1.0000 

r\ 

1.3389 

63.4350 

0.5988 

2 

1.3389 

116.5650 

-0.5988 


1.4042 

49.1066 

0.3538 

3 

1.4565 

90.0000 

-0.4855 


1.4042 

130.8934 

0.3538 


1.4014 

40.0881 

0.2256 

A 

1.4605 

73.4273 

-0.3410 

4 

1.4605 

106.5727 

0.3410 


1.4014 

139.9119 

-0.2256 


1.3810 

33.8783 

0.1540 


1.4369 

62.0404 

-0.2538 

5 

1.4521 

90.0000 

0.2904 


1.4369 

117.9596 

-0.2538 


1.3810 

146.1217 

0.1540 


Notice that if equation (9- 10a), with 0 = 0 P , is divided by (9- 14), the result is 


(9-22) 


K 1 

R n p n + l' 


Thus, if a saddle point exists, then a point of absolute minimum ( h = 0) always exists between 
the saddle point and the origin. Furthermore, Ro = 0 when P,J (cos0)= 0 and for n > 1, the values 

of 0 where this occurs (0 = 0 O ,„) are given by 0 o , n = 0p,(«-i) , where 0 p ,(„-i) signifies the Q p for n - 1 
in Table II. Since h 2 as R 0 unless R = Ro, and since the only saddle points of h 2 appear 
at the ( R P ,9 P ) in Table II, for each n, then there is no way for a particle to get from the minimum 
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associated with one value of d p to another, unless it first go out through one saddle point and then 
in through another. 

To further illustrate these ideas and as an example of the use of Table II, consider n = 3. If y > 0, 
saddle points exist at Q p = 49° and 131°. If 0 < y< 0.3538, the particle can pass through the 
saddle point at 0^, = 49° from a region where R < R p to the unbounded region where R > R p ; it 
could then enter through saddle point at 0^ = 131° into another region where R< R p . However, if 
y > 0.3538 and R < R p , the particle cannot pass through either saddle point and is trapped. 
(Remember that the particle can only move in those regions where h 2 < 1, so that these are 
allowed regions, while those for which h 2 1 are forbidden regions. Also, once a particle leaves 
a saddle point, it tends to get farther and farther away from the point multipole.) 

In order to plot h 2 we choose cylindrical coordinates; the conversion of h, given by (9-9), from 
spherical polar to cylindrical coordinates is given in Appendix B. Using these results, leads to 
Figures 1, 2, 3, 4, and 5, which illustrate allowed and forbidden regions for dipole, quadrupole, 
octupole fields, hecadecupole, and fifth-order pole for various values of y. The values of y chosen 
for Figure 1 correspond to the values used by Stormed (in his Fig. 123). For an object in space 
possessing a magnetic field well approximated by a point multipole, particles in the local space 
environment can be trapped if the associated value of y is outside the interval (y c ', y c + ), where y c ~ 
is the largest negative y and y c + is the smallest positive y for each n in Table II. 


Actual trajectories can also be numerically produced and compared to the global solutions shown 
in Figures 1 through 5. We again choose to use cylindrical coordinates, since these are easier to 
plot. We take equations (8-1) to (8-3) and use (9-8) and (C-l) to turn these equations into the 
dimensionless form [where s = vt and R = (z 2 + p 2 ) y ]: 


(9-23) 

(9-24) 


d2p -h 

(n + l)y 

3_ 

r/>‘(z/r i )j 

ds 2 ' • 

L p 2 

3p 

R n + l 

\ M 


d 2 z 3 f PjfzR- 1 )) 

ds 2 n a z r +l 

V / 


(9-25) 12!:=: 


d<p _ h n 
ds p 


The explicit form of the functions h„ is given in the table in Appendix 
In the case of a dipole field, n = \, and equations (9-23) to (9-25) take 


(9-26) 


O - h 

ds 2 ~ h{ 


2y_(2p*-z 2 ) 


R 5 


B. 

the form 
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(9-27) 


d 2 p 3pz 


(9-28) 


d§ _h i 2y 1 

ds p p 2 R 3 


These can be integrated to produce particle orbits for any value of y. Here, this was done with a 
3 rd Adams-Bashforth method 13 . The results of doing so are shown in Figure 6 for y= 1.016, and 
in Figure 7 for y = 0.999. Allowed and forbidden regions are shown in Figure 1 for y = 1.016, 
and in Figure 8 for y = 0.999. In the case where y = 1 .01 6, there are two disjoint allowed regions, 
so that a particle may be in a trapped orbit, while in the case y = 0.999, there is only one 
connected allowed region. 


10. Conclusion 


In this work, a number of novel developments have been presented. First, the magnetic vector 
potential equivalent to a magnetic scalar potential was derived. This allowed the use of 
Lagrangian dynamical methods, which require a magnetic vector potential in their formulation. 
Next, the global structure the allowed and forbidden regions of motion for an arbitrary 
axisymmetric multipole of order n was found. Explicit results were found up to order n = 5, 
while higher order results can be found by applying our general theory to n > 5, as desired. 
Finally, the topologies of the allowed and forbidden regions for n up to 5 were shown in 
graphical form. 

These results may be of use in describing the radiation belt structure of astrophysical objects that 
have insignificant magnetic dipole fields but instead have some higher order axisymmetric 
multipole structure. In the design of plasma propulsion systems that utilize magnetically 
confined plasma, certain design requirements may constrain the net magnetic field to be an 
axisymmetric quadrupole (or higher order) field, rather than a dipole field (This constraint may 
arise due to the unacceptable levels of torque that occurs when a magnetic dipole is placed in a 
planetary magnetic field). 

In the future, we plan to extend the analysis presented here for axisymmetric multipoles to the 
case of axisymmetric collections of current loops. This will allow a mapping of allowed and 
forbidden regions of motion in regions interior to a plasma propulsion device, for example. The 
magnetic field associated with a current loop is mathematically described in terms of elliptical 
integrals and functions derived from them, and the associated analysis is more difficult than the 
analysis associated with spherical harmonics. However, analysis of the motion interior to a set of 
current loops is an important one and the plan is to treat this problem as priorities, resources, and 
time allow. 
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Appendix A 


Table A Associated Legendre Functions 

Some associated Legendre functions P (x) , with x = cos0 and s = sin0 = (1 - x 2 ) 1 ' 2 , are 


n 

m 

CM 

0 

0 

1 

1 

0 

X 

1 

1 

- s 

2 

0 

(3x 2 - l)/2 

2 

1 

-3 xs 

2 

2 

3s 1 

3 

0 

x(5x 2 - 3)/2 

3 

1 

-3s(5x 2 - l)/2 

3 

2 

15x5 2 

4 

0 

(35x 4 - 3 Ox 2 + 3)/8 

4 

1 

-5xs{lx 2 - 3)/2 

4 

2 

155 2 (7x 2 -l)/2 

5 

0 

x(63x 4 - 7Qx 2 +15)/8 

5 

1 

-15^(21x 4 -14x 2 + l)/8 

6 

0 

(23 lx 6 - 3 1 5x 4 + 1 05x 2 -5)/l 6 

6 

1 

-21x5(33x 4 -3Qxr 2 + 5)/8 
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Appendix B 


The function h of the multipole pseudo-potential Jr_ in cylindrical coordinates 


„ . /n , (n + l)y P. (cos9) 

Equation (9-9): h = - + — ; — 

tfsinG R n+l 

Coordinate transformation (z, p, R are measured in terms of r s ): 

(B-l) a) z = RcosQ b)p=i?sin0 c) R = (z 2 + p 2 ) ' A 

Table B . The multipole function h = h n in cylindrical coordinates 


n 

hn 

i 

2 Y__P_ 
p i? 3 

9 


3y 3pz 

z 


p fl 5 

9 


4y 3p(4z 2 -p 2 

J 


p 2 R 1 

4 

5y 5pz(4z 2 -3p 2 ) 
p 2 tf 9 

5 

6y 

15 p(8z 4 -12p 2 z 2 + p 4 ) 

P 

8 /?" 
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Figure 1. Allowed (white) and forbidden (gray) regions for an axisymmetric 2 nd order pole 
(dipole); the polar axis is aligned with and centered on the left vertical edge of each figure. 
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T * 0 6200 


y = 0.5988 




Figure 2. Allowed (light) and forbidden (dark) regions for an axisymmetric 2 nd order pole 
(quadrupole); the polar axis is aligned with and centered on the lower horizontal edge of each 
figure. The vertical axis is p and horizontal one is z. For -y, figures reflect about z = 0. 
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y— 0.353# 


y= 0.2000 


T — 0.1000 


y - -0.4500 


y= 0.5000 


Figure 3. Allowed (white) and forbidden (gray) regions for an axisymmetric 3 rd order pole 
(octupole); the polar axis is aligned with and centered on the lower horizontal edge of each 
figure. 
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Y - 0.3600 


Y - 0.34V9 





Y - U. 1 ouo 


T ~ O.OOOU 



0.0 


-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



2.0 
15 
1.0 
0.5 
0 0 


•1.5 1.0 *0.5 0.0 0.5 1.0 I S 


Figure 4. Allowed (light) and forbidden (dark) regions for axisymmetric 4 th order pole (hexa- 
decupole); polar axis is aligned with and centered on the lower horizontal edge of each figure. 
The vertical axis is p and horizontal one is z. For -y, figures reflect about z = 0. 
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T - 0.3000 y = 0.2900 



- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 - 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 


Y = 0.2500 



- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 


2.0 
1.5 
1.0 
0.5 
0.0 

- 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 


Y - 0.1540 



Y = 0.1500 Y = t> - 0000 



- 1.5 - 1.0 - 0.5 0.0 0 5 1.0 1.5 - 1.5 - 1.0 - 0.5 0.0 0.5 1.0 1.5 


Y = -0.2250 



• 1.5 - 1.0 - 0.5 0.0 0.5 


Y = -0.2540 



0.0 


- 1.5 - 1.0 - 0.5 0.0 0.5 10 1.5 


Figure 5. Allowed (light) and forbidden (dark) regions for axisymmetric 5 th order pole; the polar 
axis is aligned with and centered on the lower horizontal edge of each figure. 
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Dipole field: y = 1.106 


1.0 i 


0.9 - 



0.1 - 


i 1 1 0.0 r 1 1 1 

-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 

z 


Figure 6. Particle motion in a dipole field for y= 1.016. 
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Dipole field: y = 1.016 



Figure 7. Particle motion in x-y plane for a dipole field with y = 1 .016. 
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• 4 


Dipole field: y = 0.999 



z 

Figure 8. Particle motion in a dipole field for y= 0.999. 



- 0.3 - 0.2 - 0.1 0.0 0.1 0.2 0.3 


z 

Figure 9. Dipole, allowed (light) and forbidden (dark) regions for y = 0.999. 
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